Thermalization of Color Gauge Fields 
in High Energy Heavy Ion Collisions 



Aiichi Iwazaki 

International Politics and Economics, Nishogakusha University, 
Ohi Kashiwa Chiba 277-8585, Japan. 
(Dated: Dec. 10, 2007) 

We discuss quantum mechanical decay ol color magnetic field generated initially at high energy 
heavy ion collisions. It is caused by Nielsen-Olesen unstable modes and is accomplished possibly in 
a period < Ifm/c. We show that the decay products, i.e. incoherent gluons may be thermalized in a 

■ sufficiently short period (< Ifm/c ). The precise determination of the period is made by caiculating 
' two point function of the color magnetic field in a model of color glass condensate. 

o : 

. In high energy colHsions of heavy ions in RHIC or LHC, the most important ingredients for producing dense quark 
^ gluon plasma ( QGP ) are smaU x gluons in the nuclei .1^, 2, 3, 4, 5]. The gluons with small transverse momenta are 
sufficiently dense in the nuclei and so they may be treated as classical fields produced by glassy large x gluons even 
after the collisions of the nuclei. It has been shown that longitudinal color magnetic and electric fields of the small 
' X gluons are generated initially at the collisions. They are classical fields and evolve classically according to a model 

■ of color glass condensate (CGC)@. It is expected phenomenologicallyQ that the decay of the gauge fields within a 
period of Ifm/c after the collisions leads to thermalized QGP. 

In an extremely high energy collision, the radius of nuclei is nearly zero due to the Lorentz contraction. Thus, 
O-i the initial gauge fields have only transverse momentum perpendicular to the collision axis, but have no longitudinal 
yi^ ^ momentum ( rapidity ). Such classical gauge fields can not possess any longitudinal momentum in their classical 
^ ' evolution, since equations of motion of the gauge fields are invariant under the Lorentz boost along the collision 
f~i , axis[l, H, y]- It has recently been shown fsl that the addition of small fluctuations with rapidity to the initial gauge 
field induces exponentially increasing modes with longitudinal momentum. Such fluctuations may be produced by 
quantum fluctuations around the classical gauge fields]y, [l3| • The production of the exponentially increasing modes 
^ implies that a process toward thermalization has started; the decay of the gauge field and the isotropization of 
momenta. Although the modes increase exponentially in time, the time needed for them to grow sufficiently large is 
, too long to be consistent with the phenomenological expectation. 

' In this letter we show that the decay of the color magnetic field is caused quantum mechanically by Nielsen-Olesen 
" [ , unstable modesflll, [l3|- It can occur in a sufficiently short period after the collisions so that the phenomenological 
OsJ ' expectation is satisfied. The decay of the gauge field produces much dense incoherent gluons. They carry longitudinal 
momenta as well as transverse momenta; their typical momenta are given by the square root, a/ gB, of a color 
magnetic field, gB, averaged over transverse plane. It is shown with reasonable choice of parameters that such dense 
gluons have sufficiently short mean free path to be thermalized immediately after their productions. Consequently, 
we obtain isotropic thermalized QGP. The life time of the gauge field is given in terms of a saturation momentum 
used in a model of CGC, e.g. the Mcleran-Venugopalan ( MV ) modelflj. It can be obtained by the calculation of 
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5^ ' a two point function of color magnetic fields in arbitral proper time, r = yjt"^ — x\ 

I r = 0. Here, without explicitly performing the calculation, we simply give general formulae for the quantities such 
as number density of produced gluons, their mean free path and temperature of QGP. All of them are expressed 
in terms of a single parameter, QgTc, Qs is saturation momentum used in the MV model and is the quantum 
mechanical life time of the magnetic field. Thus, once we obtain the value of QsTc, we can determine whether or 
not the thermaliztion of QGP is achieved at r = Tc. Indeed, when we take phenomenologically reasonable values, 
Qs = 2GeV and Tc — 0.5fm/c, or QgTc > 5, thermalized QGP is realized at t = Tq- In this paper we only discuss 
the fate of the color magnetic field and assume that the color electric field is screened immediately by dense small x 
gluons. 

First, we sketch briefly Nielsen-Olesen unstable modes generated under the color magnetic field and the subsequent 
decay of the field using SU(2) gauge theory. The modes arise in the presence of a homogeneous color magnetic field, 
Without loss of generality, we may point the field into the direction of third axis both in color space, a, and real 
space, i; Bf — BS'^'^Si^s. Then, we decompose the Lagrangian of gluons with the use of the variables, "electromagnetic 
field" Af^ — and "charged vector field" = (^4^ -I- iA'j^)/^/2 where indices 1^3 denote color components. 
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+ ig{d^A, - a,A^)$t + ^(<i>p$t - $,<i>t )2 (1) 

with = + igA^, where we have used a gauge condition, D^^^ — 0. We note that the third term in the right 
hand side represents anomalous magnetic moment of the field, <i>p. 

Using the Lagangian we find that the energy E of the charged vector field <i>p cx e^'^* in the magnetic field, 
B ~ eijdiAj, is given by E'^ = + 2gB{n + 1/2) ± 2gB. Here, ±2gB ( the integer n > ) denote contributions of 
spin components of <i>p ( Landau levels ) and denotes momentum parallel to the magnetic field. The term of i.2gB 
represents the contribution of the anomalous magnetic moments. 

Obviously, the modes with E'^{n = 0, fc| < gB) < are unstable; their amplitudes increase ( decrease ) exponentially 
in time. We call them Nielsen- Olesen unstable modes. Their presence implies that the color magnetic field is unstable; 
it decays with the production of these unstable gluons. The decay width, F, in unit volume has been calculated such 
as r = ^^^J by assuming a homogeneous color magnetic field, B. 

We can show[T3| that main results on the classical evolution of gauge fields in the Ref.Q are understood in terms of 
these Nielsen- Olesen unstable modes. Namely, the exponentially increasing modes found in the Ref. Q can be identified 
as Nielsen-Olesen unstable modes. Actually, by solving linearlized equation of motions f4j for the unstable modes 
under both homogeneous longitudinal color magnetic and electric field taken appropriately, we can find the foUowings; 
f) maximum longitudinal momentum, t'max, excited increases linearly with time, r, 2) longitudinal momentum, v, 
giving the maximum Fourier component of longitudinal pressure is small and increases very slowly with r, and 3) the 
maximum Fourier component of the pressure increases in such a way as exp(r'^/^), although it increases as exp(T^/^) 
in the Ref.Q. ( Numerical results in the Ref. Q can be fitted even by using this increasing function exp(T^/^). ) The 
saturation of the exponentially increasing component of the pressure arises due to the quartic interactions of $^ in our 
analysis. As has been recognized, the isotropization of momenta in the classical evolution needs too much time, Tjso, 
to be consistent with the phenomenological expectation, Tiso < 1 fm/c. Therefore, it is natural to see their quantum 
effects by using the above decay formula. All of the unstable modes contribute to the decay quantum mechanically. 

The unstable gluons of the exponentially increasing modes occupy the lowest Landau level ( n = ). Their wave 
functions are characterized by angular momentum, to = 0, 1, 2, , , around X3 axis, 

a;3, t; m, fcg) = (*i - i^2)VT/2 = ^^z" e^^ksxa ~ |z| + t^gB - ks) (2) 

with z ^ xi + 1x2 and normalization factor, g'^^ = ^,„i(2?^)"'+i • Here we chosen the central gauge, Ai — 
{—BX2/2, Bxi/2,0). In this formula I = 1/^/gB denotes magnetic length representing transverse extension of the 
states. The factor of exp(i-\/ gB — /c|) represents relative production ratio of the modes with each longitudinal mo- 
mentum, k^. All of the states specified by m are degenerate in energy. It has been shown [l5l.[lQ] that the states screen 
the original magnetic field in a way just like a Landau diamagnetism. Hence, the color magnetic field decays with the 
production of the Nielsen-Olesen unstable modes. This is very similar to the case that an electric field decays with 
the production of charged particles, which screen the electric field. 

It is obvious that both transverse and longitudinal momenta of the produced gluons are typically given by l/l — 
y/gB. Although the color magnetic field is homogeneous in the longitudinal direction, the gluons have the nontrivial 
longitudinal momentum. This is sharply contrasted to the case of classical evolution of the gauge fields. In the 
classical evolution, longitudinal momentum never arises dynamically due to the symmetry of initial condition and 
equation of motion of gauge fields. 

The collisions between a nuclei running into +2:3 direction and a nuclei running into —2:3 direction are supposed to 
occur at t = X3 = 0. In the high energy collisions of identical heavy ions, a longitudi nal gaug e fiel d of sm all x gluons 
is generated initially at the collisions. The explicit form of the gauge field sX t — ^Jt'^ — x\ = 2\/ x^x~ = is given 
in the MV as. 



Aj(t = 0,a:t) = ai i(a;t) -|- a2.i(xt) and as,t{xt) = -U s(xt)diUl{xt), (3) 

a 

where gauge potentials, ag^iixt) are produced by the sources, psixt^z^^ of large x gluons in s-th nuclei, l}\{xt) = 

Pexp{—ig dz~ -^pi{xt,z~)5{x~)) and U2{xt) = Pexp{—ig dz'^ -^p2{xt, z'^)S{x^)) where denotes deriva- 
tive in transverse coordinates, Xt = {xi, X2). The distribution of the sources is given such that {ps{xt, x^)ps'{yt, y^)) — 
5s,s'Ql/ g'^5'^{xt — yt)5{x^ — y^). The gauge potentials represents longitudinal color magnetic field, B^{xt) — 
igeij[ai^i(xt),a2,j(xt)\- Since this magnetic field has typically a transverse momentum, Qs, it is never homogeneous 
in transverse directions. 



3 



Although the color magnetic field is not homogeneous in space and time, it may decays due to the excitation of 
Nielsen-Olesen unstable modes mentioned above. We need to find an effective homogeneous magnetic field relevant to 
the decay in order to use the formula of the decay width, F, in the above. It is reasonable to suppose that such a field 
is given by a color magnetic field averaged over the transverse space. In order to obtain the average field explicitly, 
we use root mean square of color magnetic flux, CI{t) = (PxtB'i{xt, t). 



gB{T) = with VL^{t) = Tr(17(T)f2(T)) and (r2(r)) = 0, (4) 

where the average, {Q) is taken over the sources, ps- The integral, j^cPxt, is taken over a transverse region with 
collision area, S, relevant to a heavy ion collision, e.g. S = irRfj with nuclear radius, Rn, for central collisions. Since 
the correlation length of the magnetic field is given roughly by Qj^, we may rewrite fl^ such as D,^ = where the 

numerical coefficient c(r) can be obtained by numerical calculation of the correlation function, Tr(i?3(a;t, r)-B3(j/t, r)) 
in the formula. Thus, the homogeneous average magnetic field can be rewritten as gB — Qs\/c{t)/S. Since the 
correlation function, Tr(i?3(a;f , r)_B3(xf , r)), has been shownQ to decrease as l/r for QsT > 1, we may speculate the 
similar decrease of Tr(i?3(xf , r) 53(1/4, t)) ~ l/r. It means that c(t) ~ 1/QsT. ( It is easy to see that perturbative 
solutions, B^{xt,T), of gauge fields behaves such as B^^xt^r) oc 1/-\/t. Thus, TT{B3{xt,T)B3{yt,T)) ~ l/r for Q^t > 1. 
) Note that since c(r) does not involve a parameter with mass dimension except for Qs, it depends on t only through 
the combination, QsT. 

We assume this field as an effective homogeneous color magnetic field generating the unstable Nielsen-Olesen modes. 
Actual decay of the inhomogeneous magnetic field is more complicated than the simple case discussed in the present 
paper. But, by using the effective magnetic field we may obtain approximately the life time of the field and the 
number density of gluons produced subsequently. In order to determine more accurately these quantities, we need to 
examine more closely how the inhomogeneous color magnetic field decays. 

We should make a comment on a consistency condition that Nielsen-Olesen unstable modes have to be located inside 
of the region involving the magnetic field. This constraints the parameter, gB, such as Xj \fgB < ^ S/ir. Namely, 
the magnetic length, 1/ \/ gB, is smaller than the radius, ^ S/tt, of the collision area. This condition is satisfied in 
general even for peripheral coUisions as far as QlS > tt'^ /c{t) with Qs> I GeV. 

As we can see soon below, the decay of the magnetic field is very rapid. On the contrary, the dependence of gB on 
T is much smooth ; gB = Qs\/c{t)/S oc r~^/^. Hence, we may apply the above decay formula to gB depending on 
T. Consequently, we find the life time of the magnetic field by solving the following equation, 



i = r.i?... = M)!M = £(l^, (5) 

where — t denotes longitudinal extension of the color magnetic field. ( In the explanation of Nielsen-Olesen 
unstable modes we have used standard time coordinate, not r = ^Jt'^ ~ Thus, we can apply the decay formula 
of B only to midrapidity, 77 — log(^/^|^^) ~ 0. ) It is important to note that does not depend on transverse area, 

S. It implies that the life time of magnetic field does not depend on the centrality of the collisions. 

Solving the equation®, we obtain the value of QsTc- Thus, if Qs becomes larger, the life time, r^, of the color 
magnetic field becomes smaller. Since the saturation momentum, Qs, is larger in LHC than that in RHIC, QGP is 
expected to be thermalized much earlier in LHC that that in RHIC. 

Although the precise evaluation of c(t) is needed for solving the equation, we may estimate roughly Tc by assuming 
a naive correlation, lT:{B^[xt)Bz{yt)) = Q^jg^ for \xt — yt\ < Qj^ and vanishes for otherwise. Then, it leads to 
il^ ~ Q^S/g^, namely, c = 1. Hence, we find Tc ~ 5Q7^ ~ 0.5 fm/c for Qs = 2GeV. This satisfies the phenomenological 
requirement, Tc < Ifm/c. Although this estimation is very rough, the precise evaluation of c{tc) would similarly give 
a phenomenologically reasonable value in the following way. 

Since the real correlation, TT{B3(xt) B^{yt)) , has no such sharp boundary like exp(— Qgr) in r = \xt —yt\ as assumed 
above, the integral J (Pxt gives much larger contribution than one estimated just above. This results in c(t = 0) 3> 1. 
On the other hand c(t) decays as c(t) oc 1/QsT. These would lead to c(tc) ~ 0{1). Thus, we expect that an 
appropriate result, e.g. Tc = 0.5 fm/c, would be also obtained even with the precise evaluation of TT{Bz(xt) B3{yt)) . 

We now proceed to show that gluons produced by the decay of the magnetic field are sufficiently dense to be 
thermalized quickly when a phenomenological requirement, Tc < 1 fm/c, is satisfied with 1 GeV < < 3 GeV. 
We give a general formula for the density of gluons and determine whether or not the gluons are thermalized. The 
thermalization is assumed to be achieved when a produced gluon interacts with others many times within a period. 
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Tc. We also give their temperature when the gluons are thermahzed. The formulae are expressed in terms of the 
parameter, QgTc- 

The momentum distribution of the gluons produced by the decay is proportional to absolute square of the wave 
function of the unstable gluons, that is, |<&(A;t, /ca, t; m)p = | J (Pxt exp(— z/ctXt)$(a;t, 2:3, t; to, /ca)^. 
Hence, the gluon number density is given by 

Nir) ^ No{t) f d^ktdks \Hkuh,T;m)\^ (6) 

•' 0<m<NR 

with normalization constant, Nq(t), where we have summed equally over the wave functions with different m since 
they are degenerate. The largest angular momentum, Nn is determined by the condition that the gluons must be 
located inside the transverse region with surface area of S: Because the average ($||zp|<i>) — 2{m + 1)P should be 
less than S/tt, Nr = ^{R/l f - 1 ~ \{R/lY = V^SQ^/QsT with ttR"^ = S. On the other hand, the normalization 
constant, No{tc), is determined in the following. Namely, the energy density of the produced gluons is equal to the 
energy density of the color magnetic field, B^{t)/2 = Tr((B(xt, r)B(xt, r)))/2. 



= 7Vo(Te) / d^hdh^k^ + ki Mkt,k3,T,;m)\'', (7) 

" ^ AT - 



^ •' 0<m<A'R 



where we have used the energy, ^Jk^ + A:| , of free gluons because the gluons become almost free after the decay of 
the magnetic field. We have also assumed the translational invariance in the transverse directions. Thus, Nq(tc) is 
given by 

N,ir.)^ (£^MZ . (8) 

252 / d^hdksy/k^ + klY.o<m<Nn Mh, k^, T,; m)\^ 

Inserting the formula into eq®, we find the number density of the gluons, 

, (.gB(re))2 Jd^ktdksZ 0<rn<NR\^ikt,k3,Tc;m)\'^ 

~ V / d^hdk^y/k-^ + k^ Eo<m<N,, Mku fc3, T,; to)P 

Noting that 

mkt,k3,T;m)\^ = \AfTrg^{2if)"'{ki + ik^T (^M^k^x^ ~ k^P + T^gB - kl)f (10) 

= SPtt ^^^'?" exp(-2/2fc,2 + 2r\lgB - k\) (11) 
to! V 

= SZ^TT^^^^ exp(-2(72 + 2fy^r^) (12) 
with l'^ — l/gBir), at = Ikt, era = Zfcs and r = t/1, we obtain 

^ 1 {gB{r^)r /|,3|<,rf2a.da3E:^to^exp(^2a,2 + 2f.yr^) 

dir.)Ql^.fSQl\^'^^^^^^^^^^ (14) 



2^2 \^ Stt 

with igB{Tc)f = d{Tc)Qt and = {8tt/SQI)^^^^/Q^, where 



■ exp(— 2o-2 -I- 2tc\/1 - (T; 



L3K1 d^a.daav/^FT^ E™io ^ exp(-2a2 + 2f,yr^) 



(15) 
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Here we note that the value of d(r) = {gB{T))^ /Q'^ has already been calculated in the referencefs'l such as ^(t) ~ 
0.15(1/Q,r) for Q,t > 1. 

Up to now, we give a general formula of the gluon density in terms of a quantity, QsTc or c{tc) — Sn / (Q sTc)'^ , which 
can be obtained by calculating the quantity, J (fxtd^ytTT{B3{xt,T)B3{yt,T)) = Hill^i^ in the MV model. We should 

note that the parameters in N{tc), such as c{tc) = %tt / {Q sTc)'^ , Nr = yj2SQ1/QsTc and fc = (STr/SQ^)^/^ VQsT'c, 
are all expressed in terms of QsTq- Thus, once we determine the parameter, QsTc as well as Qs and g^, we can obtain 
the number density of the gluons produced aX t — Tc- 

Although we have not yet determined QsTc by calculating explicitly c(r), we will estimate how large is the number 
density, N{tc) and determine whether or not the gluons are thermalized. Assuming a phenomenologically expected 
value of Tc = 0.5 fm/c with Qs = 2 GeV and R = 6.4 fm, we obtain N{tc) ^ 21/fm'^ with g — 2. Since the 
number density is very large, the gluons may be thermalized quickly. Actually, the cross section, a, of the gluons with 
the typical momentum, \/gB ( = Qs{VQsTc)~^{8tt/SQI)^^'^ = 188 MeV), is approximately given by cr = ag/gB ( 
Qfg = g^/47r ). Hence, the mean free path, l^, is that /,„ = 1/(jN{tc) ~ 0.06 fm, where we have used ag — 1 since 
the typical momentum of the gluons is y/gB = 188 McV. Thus, the number of collisions for a gluon to interact with 
others within the period, Tc = 0.5 fm/c, is given such that Tcc/lm — 8 times. This may be sufficiently large that the 
gluons are thermalized within the period, Tc — 0.5 fm/c. Therefore, after the decay of the color magnetic field, much 
dense gluons are produced and are thermalized immediately. 

We have assumed the parameter such as QsTc = 5 in the above estimation. If we assume a larger one, e.g. 
QsTc = 10 with Qs fixed, the number of collisions within a period, Tc becomes larger and so the thermalization 
proceeds more efficiently. This is because the number of collisions behaves as Tcc/lm (x ^JQs^c since N{tc) cx I/^/QsTq 
and a cx QsTq- While, when we increase the saturation momentum, Qs with QsTc fixed, the number of collisions 
increases as Tcc/Zm c>c qV^ since N{tc) (x qI^^ and a oc Qj^. Hence, the thermalization proceeds more efficiently in 
higher energy collisions of heavy ions. Similarly, when we take larger collision area, SQ^ than the previous case, the 
thermalization proceeds more efficiently since Tcc/lm (oc (SQl)"^^^). This implies that when the collision is central, 
the most efficient thermalization is achieved. ( The dependence of f{Nfi,fc) on QsTc is not so serious to change 
significantly the above results. ) 

Finally, by taking the parameters used above such as = 2 GeV, Tc — 0.5 fm/c, we estimate the temperature, Tc, 
of the QGP by assuming a gas of free quarks and gluons, 

(-^ + "f^^^c " ^2 d{Tc)Qs/2g (16) 

with gluon's ( quark's ) number of degree of freedom, rig = 16 ( = 36 ). Since d{QsTc = 5) ~ 0.03, we find that 
Tc ~ 246 MeV. 

To summarize, we have shown that the coherent color magnetic field generated initially in high energy heavy ion 
collisions decays to produce incoherent dense gluons. The decay occurs owing to the quantum effects of Nielsen- 
Olesen unstable modes. The life time of the magnetic field has been given in terms of the quantity, c(t), which is 
obtained by the evaluation of Ti{B3{xt, T)B3{yt, r)) in the MV model. We have expressed the number density of the 
produced gluons and their mean free path in terms of QsTc Assuming Tc = 0.5 fm/c with the saturation momentum, 
Qs ~ 2 GeV, we have shown that the produced gluons can be thermalized within the period, Tc. When the saturation 
momentum becomes larger, the thermalization is achieved earlier. In our discussion inhomogeneous color magnetic 
field has been assumed to decay with the production of the Nielsen-Olesen unstable modes. They are supposed to be 
induced by an effective homogeneous color magnetic field given by the root mean square of the original color magnetic 
flux. In order to see more precisely the decay, we need to examine more closely how the inhomogeneous color magnetic 
field decays generating the unstable modes. 

We would like to express thanks to Dr. Itakura for useful comments. 



[1] L.D. McLerran and R. Venugopalan, Phys. Rev. D49, 2233 (1994); D49, 3352 (1994); D50, 2225 (1994). 

[2] A. Krasnitz and R. Venugopalan, Nucl. Phys. B557, 237 (1999); Phys. Rev. Lett. 84, 4309 (2000); 86, 1717 (2001). 

[3] A. Krasnitz, Y. Nara and R. Venugopalan, Phys. Rev. Lett. 87, 192302 (2001); Nucl. Phys. A717, 268 (2003). 

[4] T. Lappi, Phys. Rev. C67, 054903 (2003); C70, 054905. 

[5] T. Lappi and L. McLerran, Nucl. Phys. A772, 200 (2006). 



6 

[6] E. lancu, A. Leonidov and L. McLerran, hep-ph/0202270' 

E. lancu and R. Venugopalan, hep-ph/0303204 
[7] T. Hirano and Y. Nara, Nucl. Phys. A743, 305 (2004); J. Phys. G30, S1139 (2004). 

[8] P. Romatschke and R. Venugopalan, Phys. Rev. Lett. 96, 062302 (2006); Phys. Rev. D74, 045011 (2006). 
[9] K. Fukushima, F. Gelis and L. McLerran, Nucl. Phys. A786, 107 (2006). 
[10] F. Gelis, T. Lappi and R. Venugopalan, hep-ph/0708.0047. 

[11] N.K. Nielsen and P. Olesen, Nucl. Phys. B 144 (1978) 376; Phys. Lett. B 79 (1978) 304. 
[12] G.K. Savvidy, Phys. Lett. B 71 (1977) 133. 

H. Pagels, Lecture at Coral Gables, Florida, 1978. 
[13] Hokkaido meeting. The Physical Society of Japan, 2007, in preparation. 
[14] K. Fukushima, Phys. Rev. C76, 021902(R) (2007). 
[15] H.B. Nielsen and M. Ninomiya, Nucl. Phys. B156, 1 (1979). 

H. B. Nielsen and P. Olesen, Nucl. Phys. B160, 330 (1979). 
[16] A. Iwazaki, hep-ph/0705085^ 



